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Abstract. We express a weighted generalization of the Delannoy numbers in terms of 
shifted Jacobi polynomials. A specialization of our formulas extends a relation between 
the central Delannoy numbers and Legendre polynomials, observed over 50 years ago, to 
all Delannoy numbers and certain Jacobi polynomials. Another specialization provides 
a weighted lattice path enumeration model for shifted Jacobi polynomials, we use this 
to present a new combinatorial proof of the orthogonality of Jacobi polynomials with 
natural number parameters. The proof relates the orthogonality of these polynomials 
to the orthogonality of (generalized) Lagucrre polynomials, as they arise in the theory 
of rook polynomials. We provide a combinatorial proof for the orthogonality of certain 
Romanovski- Jacobi polynomials with zero first parameter and negative integer second 
parameter, considered as an initial segment of the list of similarly transformed Jacobi 
polynomials with the same parameters. We observe that for an odd second parameter one 
more polynomial may be added to this finite orthogonal polynomial sequence than what 
was predicted by a classical result of Askey and Romanovski. The remaining transformed 
Jacobi polynomials in the sequence are either equal to the already listed Romanovski- 
Jacobi polynomials or monomial multiples of similarly transformed Jacobi polynomials 
with a positive second parameter. Wc provide expressions for an analogous weighted 
generalization of the Schroder numbers in terms of the Jacobi polynomials, and use this 
model, together with a result of Mansour and Sun. to express the Narayana polynomials 
in terms of shifted Jacobi polynomials. 



Introduction 

It has been noted more than fifty years ago [HI EHl El] , that the diagonal entries of the 
Delannoy array {dm.n), introduced by Henri Delannoy [7], and the Legendre polynomials 
Pn{x) satisfy the equality 

(1) dn,n = Pn(3), 

but this relation was mostly considered a "coincidence". The main result of our present 
work is that a weighted lattice path enumeration model may be used to show that ([I]) can 
be extended to 

(2) dn+a,n = Pt'°\^) for all « G Z such that «> -n, 
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where Pn ' (x) is the Jacobi polynomial with parameters (a,0). This result indicates 
that the interaction between Jacobi polynomials (generalizing Legendre polynomials) and 
the Delannoy numbers is more than a mere coincidence. We then use this model to 
provide a new combinatorial proof of the orthogonality of the Jacobi polynomials with 
nonnegative integer parameters. The first combinatorial proof of this orthogonality is 
due to Foata and Zeilberger [TU], who provided a combinatorial model for the weighted 
integral of an arbitrary product of Jacobi polynomials with arbitrary parameters a, (3. The 
weighted lattice path enumeration model presented here is used to prove orthogonality 
only, and we need to assume G N. In exchange, our model is "more combinatorial" 
in the sense that the parameters G N also "count something" and are not just 
"mere weights" . Using our model we may also consider Jacobi polynomials with negative 
integer parameters and show that, for a positive integer the first (3 — 1 entries in 
the sequence of Jacobi polynomials {Pn~^\x)}n>o form a symmetric sequence, while 
the entries of degree at least j3 are monomial multiples of Jacobi polynomials of the 
form P^l^i]{x). The first — l)/2j entries in this sequence form a finite sequence of 
orthogonal polynomials, discovered by Romanovski [26], although he and Askey [2] made 
the statement in greater generality but about the first [(/? — 2)/2j entries only. We observe 
that in this special case, for odd (3 the sequence may be extended by an additional entry, 
and provide a combinatorial proof of this amended statement. Finally, by extending our 
model to Schroder paths, we show that ordinary Schroder numbers too may be obtained 
by substituting 3 into an appropriate Jacobi polynomial, and prove a formula for iterated 
antiderivatives of the shifted Legendre polynomials. This formula, together with a result 
of Mansour and Sun [23], allows us to express the Narayana polynomials in terms of 
shifted Jacobi polynomials. 

Our paper is structured as follows. In Section [2] we introduce weighted Delannoy paths 
associating weight u to each east step, weight v to each north step and weight w to each 
northeast step. This model was first considered by Fray and Roselle [TT]. Using their 
formulas, we show that the total weight of all Delannoy paths in a rectangle may be 
expressed by substitution into a shifted Jacobi polynomial with the appropriate parame- 
ters. The shifting is achieved by replacing x with 2a; — 1 in the definition of the Jacobi 
polynomial. 

A specific substitution into the parameters m, f , w yields the shifted Jacobi polynomials 
themselves. Using this observation, in Section [3] we present our lattice path enumeration 
proof for the orthogonality of the Jacobi polynomials P!f'^\x) when a, (3 G N. The 
proof reduces the weight enumeration to a known formula in the classical theory of rook 
polynomials, used to show the orthogonality of the (generalized) Laguerre polynomials. 
Further discussion of this connection may be found in Section HJ 

The discussion of the properties of the sequence of Jacobi polynomials {Pn' ^■'(x)}„>o 
is in Section [5|, while the results associated to our weighted generalization of Schroder 
numbers may be found in in Section [61 The paper inspires as many questions as it 
answers, some of which are listed in the concluding Section [3 
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1.1. Delannoy numbers. The Delannoy array {dij : i,j E Z) was introduced by Henri 
Delannoy [7i|. It may be defined by the recursion formula 

(3) dij = di-ij + dij-i + di-ij-i 

with the conditions d^^Q = 1 and dij = if i < or j < 0. The historic significance of these 
numbers is explained in the paper "Why Delannoy numbers?" [3] by Banderier and Schwer. 
The diagonal elements {dn,n : n > 0) in this array are the central Delannoy numbers 
(A001850 of Sloane [22] )• These numbers are known through the books of Comtet [6] and 
Stanley [29|, but it is Sulanke's paper [33] that gives the most complete enumeration of 
all known uses of the central Delannoy numbers. 



1.2. Jacobi, Legendre and Romanovski polynomials. The n-th Jacobi polynomial 
P^'^\x) of type is defined as 

Jn 

p("'(^\x) = (-2)-"(n!)-i(l - x)-°(l + ((1 - a;)"+"(l + . 



If a and (3 are real numbers satisfying a,f3 > —1, then the polynomials {Pn°''^\x)}n 
form an orthogonal basis with respect to the inner product 

(/, g) ■■= t /(^) ■ ■ (1 - x)-(l + xf dx. 



>0 



The following formula, stated only slightly differently in [32l (4.21.2)], may be used to 
extend the definition of P^'^\x) to arbitrary complex values of a and (3: 

The following relation provides a way to "swap" the parameters a and /3 (see [S], Ch. V, 
(2.8)]). 

(5) (-l)"Pi"''3)(-x) = '")(a;). 



The Legendre polynomials {Pn(a^)}n>o are the Jacobi polynomials {Pn'^\x)}n>o- They 
form an orthogonal basis with respect to the inner product 

{f,9)-=J f{x)-g{x)dx. 
Substituting a = /? = into (jlj) yields 

As a specialization of ([5]) we obtain 

(7) (-1)"P„(-X) = Pn{x 
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The shifted Legendre polynomials Pn{x) are defined by the linear substitution 

Pn{x) :=P„(2x-l). 
They form an orthogonal basis with respect to the inner product 

{f,9) ■■= / f{x) - gix) dx. 



xK 



They may be calculated using the following formula: 

(s) ^»M^E(-r'(:)(":>'^i:(-ri;:!:)(T 

k — k — 

A generalization of ([8]) is stated and shown in Proposition [221 Shifted Legendre polynomi- 
als are widely used, they even have a table entry in the venerable opus of Abramowitz and 
Stegun [H 22.2.11]. Their obvious generalization, the shifted Jacobi polynomials Pn°''^\x), 
defined by the formula 

(9) P("'^)(x):=i'i-'^)(2x-l), 

seem to receive rare mention by this name, a sample reference is the work of Gat- 
teschi [12] . On the other hand, a finite sequence of orthogonal polynomials, discovered by 
Romanovski [221 P- 1025] is essentially the same as a slightly differently shifted sequence 
of Jacobi polynomials. Romanovski's result was rephrased by Askey [H (1.15) and (1.16)] 
as follows. 

Theorem 1.1 ( Askey- Romanovski) . When a > — 1 and a + f3 + N + 1 < then for 
< m,n < N/2 the polynomials 

= + ^)% Fi(-n, n + a + P + l;a + l;x) 
satisfy the orthogonality relation 

(_l)n+ir(_^ _ ^ _ p^Y{n + a + l)(/5 + 1)„ . S^^^ 



Rl^'f^\x)Rl^^'^\x)x"{l+x)'^ dx 



{2n + a + p + l)n\T{-p) 



Here and in the rest of this paper, Sm,n is the Kronecker delta, (a -|- 1)„ = (a -|- l)(a -|- 
2) ■ ■ ■ [a + n) is the Pochhammer symbol and T{z) is the gamma function [U 6.1.1]. The 
the Gaussian hypergeometric function 2P1 will not used beyond the end of this subsection. 

Comparing the above definition of R^^'^^ (x) with the analogous formula for the Jacobi 
polynomials found in the database created by Koekoek and Swarttouw [T^ (1.8.1)] we 
obtain 

(10) Rl^'''\x) = Pi'''''\2x + 1). 

This relation was noted by Chen and Srivastava [H (4.7)]. 

Remark 1.2. In his very concise note, Romanovski [26] considered several finite sequences 
of orthogonal polynomials. Not all finite orthogonal polynomial sequences called "Ro- 
manovski polynomials" in the literature coincide with the polynomials R^'^\x) above. 
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To avoid confusion, Lesky [22] distinguishes between Romanovski-Bessel, Romanovski- 
Jacobi and Romanovski-Pseudo Jacobi polynomials. The polynomials Ri^'^\x) above 
are examples of Romanovski- Jacobi polynomials. 

1.3. Favard's theorem. Our main reference for the basic facts on orthogonal polynomi- 
als is Chihara's book [5]. A moment functional £ is a linear map C[x] ^ C. A sequence 
of polynomials {pn{x)}'^^Q is an orthogonal polynomial sequence with respect to £ if Pn{x) 
has degree n, C{pm{x)pn{x)) = for m 7^ n, and 7^ for all n. Such a sequence 
exists if and only if C is quasi-definite (see [5l Ch. I, Theorem 3.1], the term quasi-definite 
is introduced in [5l Ch. I, Definition 3.2]). Whenever an orthogonal polynomial sequence 
exists, each of its elements is determined up to a non-zero constant factor (see [5l Ch. I, 
Corollary of Theorem 2.2]). 

A way to verify whether a sequence of polynomials is orthogonal is Favard's theorem [5l 
Ch. I, Theorem 4.4]. This states that a sequence of monic polynomials {p„,(x)}„,>o is an 
orthogonal polynomial sequence, if and only if it satisfies the recurrence formula 

(11) Pn{x) = {X - Cn)Pn-l{x) - \nPn-2{x) n = 1, 2, 3, . . . 

where p_i(x) = 0, po{x) = 1, the numbers c„ and are constants, A„ 7^ for n > 2, 
and Ai is arbitrary (see [5], Ch. I, Theorem 4.1]). Conversely, for every sequence of monic 
polynomials defined in the above way there is a unique quasi- definite moment functional 
£ such that £(1) = Ai and {P„(x)}J^q is the monic orthogonal polynomial sequence with 
respect to £. The original proof provides only a recursive description of £ for a given 
{Pn{x)}n>o satisfying (ITTi) . Viennot [M] gave a combinatorial proof of Favard's theorem, 
upon which he has built a general combinatorial theory of orthogonal polynomials. In his 
theory, the values £(a;") are explicitly given as sums of weighted Motzkin paths. 

1.4. Central Delannoy numbers and Legendre polynomials. Equation ([1]) linking 
the central Delannoy numbers to the Legendre polynomials has been known for over 
50 years [H], [20], [23]. Until recently there was a consensus that this link is not very 
relevant. Banderier and Schwer [3] note that there is no "natural" correspondence between 
Legendre polynomials and the original lattice path enumeration problem associated to the 
Delannoy array, while Sulanke [33] states that "the definition of Legendre polynomials 
does not appear to foster any combinatorial interpretation leading to enumeration" . 

The present author made two attempts to provide a combinatorial explanation. In [TB] 
we find a generalization ([T]) to substitutions of 3 into Jacobi polynomials and an asym- 
metric variant of the Delannoy array, having the same diagonal elements. In [TB] we 
find a construction of a polytope whose face numbers are the coefficients of the powers of 
(x — 1)/2 in ([6]) and the Delannoy numbers enumerate faces having nonempty intersections 
with certain generalized orthants. Both attempts are "imperfect" in some sense: the first 
does not relate the original Delannoy array to substitutions into Jacobi polynomials, the 
second does not involve lattice path enumeration. The present work overcomes both of 
these "imperfections" . 
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1.5. Rook polynomials and orthogonal polynomials. An excellent short summary 
of the classical theory of rook polynomials is given by Gessel in [131 Section 2] . For more 
information we follow his suggestion and refer the reader to Kaplansky and Riordan [18] 
and Riordan [23 pp. 163-277]. Let 5 be a subset of [n] x [n], i.e., a board. Here [n] is 
a shorthand for {1, . . . , n}. A subset S* of 5 is called compatible if no two elements of S 
agree in either coordinate. The rook polynomial of B is defined as 

n 

rsix) := 5^(-l)Vfca;"-'= 

fc=0 

where is the number of compatible fc-subsets of B. Let C be the linear functional on 
polynomials in x defined by := n\. Then 

POO 

C{p{x)) = I e~^p{x) dx 
Jo 

and the number of permutations n of [n] x [n] such that no (z,7r(z)) belongs to B is 
£(rij(x)). 

The rook polynomial of [n] x [n] is the Laguerre polynomial 

n / \ ^ 

(12) a^)^=E(-i)'u 

k=0 ^ ^ 

In terms of the simple Laguerre polynomial as usually normalized, 

Ux) = (-l)"n!L„(x). 

It can be shown combinatorially that the Laguerre polynomials form an orthogonal basis 
with respect to the inner product induced by C. We have 

C{lm{x)ln{x)) = S,n,nn\ 

Here 6m,n is the Kronecker delta. 

2. Weighted Delannoy numbers and shifted Legendre polynomials 

A Delannoy path is a lattice path using only three kinds of steps: (0, 1) (east), (1,0) 
(north), and (1, 1) (northeast). One of the most plausible generalizations of the Delannoy 
array is the following. 

Definition 2.1. Let u, v, w be commuting variables. We define the weighted Delannoy 
numbers d^^^ as the total weight of all Delannoy paths from (0,0) to {m,n), where each 
east step (0, 1) has weight u, each north step has weight v, and each northeast step has 
weight w. The weight of a lattice path is the product of the weights of its steps. 

The usual Delannoy numbers dm,n are obtained by substituting u = v = w = 1. 
Without any reference to Delannoy numbers, a short study of the numbers d^^'^ may be 
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found in the work of Fray and Roselle [HI Section 2]. They use the notation f{n, k) := 
(fyfk^. The recursion formula ([3]) of the ordinary Delannoy numbers generahzes to 

(13) d. =u- d^Lij + V ■ dijl^ + w ■ d^L{-_^, 

see [TTl Eq. (2.1)]. The weighted Delannoy numbers also satisfy the following explicit 
formula [HI Eq. (2.2)], which was shown by Fray and Roselle using generating functions: 



fc=0 ^ / ^ / 
Using the identity 

^m\fm + n — k\ ( m + n — k \ fm + n — k\ fm + n — 2k 

kj\ m J \k,m — k,n — kj \ k J\ n — k 

Equation (fT4|) may be restated as follows. 

Proposition 2.2 (Fray- Roselle) . The weighted Delannoy numbers are given by 



k^oK k J\ n-k J 



Note that Equation f[T^ and Proposition 12.21 have the following one-line proof: a De- 
lannoy path from (0, 0) to (m, n) containing k northeast steps must contain (m — k) east 
and {n — k) north steps, and these steps may be listed in ^'^^'n-k) ^^Y^- 

Setting m = n in Proposition 12.21 yields the following formula for the weighted central 
Delannoy numbers: 

(-) «-Er;*)e:::f) 

k=0 ^ ^ ^ 

Using ([8]) and (fTSl) we may show the following formula. 

Lemma 2.3. The weighted central Delannoy numbers are linked to the shifted Legendre 
polynomials by 

cr = {-^rh (--) . 

V w J 

Proof. In a Delannoy path from (0, 0) to (ra, n) the total number of north and northeast 
steps is n. Thus we have 

d'-nT = i-wrdi'-^/-'-' 

which is equal to {—w)^dn,n^^^^' ^ since the number of east steps is the same as the 
number of north steps. By f|T5|) we have 

'2n ~ k\ f2n — 2k\ / uv\'^-~k 



r:::r (-^) (-^ 



fc=0 

and the statement follows from ([8]), after replacing k with n — k. □ 
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Substituting u = v = w = 1 into Lemma 12.31 yields 

dn,n = = (-l)"Pn,(-l) = (-1)'^ " Pn,(-3). 

Equation ([1]) now follows from the "swapping rule" ([7]). We may use the same rule to 
rewrite Lemma [2.31 as follows. 

Proposition 2.4. The weighted central Delannoy numbers and the shifted Legendre poly- 
nomials satisfy 

Proof. By definition, we have (— 1)"P„(— x) = (— 1)'^P„(— 2a: — 1). Using ([7]) we obtain 

(-l)"P„,(-2x - 1) = P„(2x + 1) = Pnix + 1). 
Thus we may rewrite ([7]) for shifted Legendre polynomials as 

(16) {-irP4-x) = p^{x+i). 

The statement now follows immediately from Lemma [2.31 □ 

Remark 2.5. Using Lemma [2.31 and Proposition 12.41 we obtain two infinite sets of weight- 
ings yielding the central Delannoy numbers dn,n as the total weight of all Delannoy paths 
from (0,0) to {n,n). Indeed, Lemma [2.31 implies 

^^n,n = forallr gM\{0}, 

whereas Proposition 12.41 yields 

<n = <„/^'-' for allr gM\{0}. 

In order to extend the validity of our formulas to non-central weighted Delannoy num- 
bers, we prove the following generalization of Equation ([8]). 



Proposition 2.6. For all a and (3 ^ C we have 

n+a 
k=0 



Proof. Using (jl]) we may write 



irp(-'^){x) = j2( 

j=0 ^ 



+ a + 3 + j\ (n^ a 
\x - irt^r'^^'Kx] = > I .\\x 

3 J\n-J. 



i=o ^ J \ ~ J/ j^^Q \ 



Since 

+ a\ f <y + j\ f n + a\ fn + a — k\ fn + a\ fn + a — k 



a + jj \ k J \ k J \a + j — k J \ k J \ n—j 
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changing the order of summation in the previous equation yields 

f^fn + a\ ^ sa+i-k (n + a + l3 + j\(n + a- k 



fc=0 ^ ^ j=0 

(n+a—k\ 



3 / \ n-j 

It should be noted that ("+"~'^) = H a < k and j < k — a. Since 

( ly^^^^J^^ _ - a - (3 - I 

the last equation may be rewritten as 

n — a — p — 1\ ( n + a — k 
n-j 



n+a / _L \ / 



fc=o ^ ^ j=0 

Using the well-known polynomial identity 

we obtain 

/—n — a — j3 — l\ f n + a — k\ f~P — k — 1 

The statement now follows from the last equation for (x — l)"Pi°''^^(x) and from 

-P-k-l\^^_^^J(3 + k + n 
n J \ n 

□ 

Corollary 2.7. In the case when a = we obtain 

k=0 ^ 



/c/ \ 77, 



Using Corollary 12.71 and Proposition 12.21 we may generalize Lemma 12.31 to the main 
result of this section. 

Theorem 2.8. The weighted Delannoy numbers and the shifted Jacobi polynomials are 
linked by the formula 

Here j3 E 1^ is any integer satisfying j3 > —n. 



Proof. By Proposition 12.21 we have 



_^ (2n + (3-k\(2n + (3- 2k\ 

k=0 ^ ^ ^ ^ 
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Replacing k with {n — k) yields 

n + P + k\l'p + 2k 
n — k / \ k 



n 

"n+/3,n - / . 

Here 



n + p + k\/'P + 2k\ _ / n + p + k \ _ /n + p + k\ /n 
n-k )\ k J~\n-k,f3 + k,kJ~\ n J \k 



thus we also have 

n 

iu,v,w \ 

k=0 

On the other hand, by Corollary 12.71 we have 



)-r'(-^)-E(:)( 



n 



which differs from d^'^'^^ only by a factor of u'^. □ 
Direct substitution oi u = v = w = 1 into Theorem 12.81 yields 

(17) rf„+^,n = {-irpi'''\-i) = 

a simpler formula may be obtained by using the "swapping rule" (jSD, which, after replac- 
ing the letter /5 with a, yields ([2]). We conclude this section with the generalization of 
Proposition 12.41 to all weighted Delannoy numbers. 

Proposition 2.9. The weighted Delannoy numbers and the shifted Jacobi polynomials 
satisfy 

Proof. In analogy to the argument seen in the proof of Proposition 12.41 we may rewrite 
(12]) for shifted Jacobi polynomials as 

(18) (-l)"Pi"'^H-x) = Pf'")(a; + l). 

The statement is an immediate consequence of this swapping rule and Theorem 12.81 □ 



3. Orthogonality of Jacobi polynomials with natural number parameters 



As a consequence of Theorem 12.81 we obtain the following lattice path representations 
of the shifted Legendre polynomials P„ and the shifted Jacobi polynomials 

Pr^(x). 

Corollary 3.1. For all (3 ^Tj satisfying [3 > —n, the shifted Jacobi polynomial Pi^''^\x) 
is d^^^pl^, i.e., the total weight of all Delannoy paths from (0, 0) to {n + l3,n), where each 
east step has weight 1, each north step has weight x, and each northeast step has weight 
-1. 
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In this section we use this Corollary to provide a combinatorial proof for all a,/? G N, 
of the fact that the Jacobi polynomials P^'^\x) form an orthogonal basis with respect 
to the inner product 

(/, g) ■■= [ fix) ■ g{x) ■ (1 - a:)"(l + x)^ dx. 



u)''u'^ du 



Using the substitution u := [x + l)/2 yields 

j fix) ■ g{x) ■ (1 - + x)^ dx = 2"+^+^ ^ /(2m - 1) ■ g{2u - 1) ■ (1 

thus the orthogonality relations we want to prove may be restated for shifted Jacobi (and 
Legendre) polynomials as follows. 

Theorem 3.2. For all a, (3 gN, the shifted Jacobi polynomials P^'^\x) form an orthog- 
onal basis with respect to the inner product 

{f,9) ■■= C f{x)-g{x)-{l-x)'-x'' dx. 



First we prove Theorem 13.21 for the case a = only. Since an inner product of any 
polynomial with itself is nonzero, it is sufficient to show the following orthogonality rela- 
tion. 

Proposition 3.3. For all m,n,j3 G N satisfying m < n we have 

»i _ 

^m+l3 . p{0,P)(x) dx = 0. 





Proof. Using the lattice path representation stated in Corollary 13.11 each Delannoy path 
from (0, 0) to {n+(3, n) containing k north steps contributes a term (— to Pn ^ {x) . 
and a term {^-ly-k^k+m+p ^m+p . p^^^l^) Since 

•1 ^ 
dx 



k + m + (3 + l 

we obtain that each Delannoy path from (0, 0) to {n + /3, n) containing k north steps 
contributes a term 

*° (n + m + /? + l)!.^^™+^-p(o./5)(a:)dx. 

Therefore ^ 

{n + m + p+l)\- [ x"*+^-P^°'^)(x) dx 







is the total weight of all pairs (L, a) where L is a Delannoy path from (0, 0) to {n + /3, n) 
and (J is a bijection {r, oi, . . . , a„+/3, bi, . . . , b^} — > {1, . . . , m + n + j3 + 1} , subject to the 
following rules: 

(i) a{r) < cr(aj) holds for all i such that there is an east step in L from {i — 1, y) to 
{i, y) for some y; 
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(ii) cr(r) < (T{bj) holds for j = 1, 2, . . . , m. 

The weight of the pair (L, a) is defined to be a function of the Delannoy path L: each 
northeast step contributes a factor of (—1), all other steps contribute a factor of 1. 

Indeed, a Delannoy path L from (0, 0) to {n + /3, n) with north steps has n — k 
northeast steps and n + P — {n — k) = (5 + k east steps. After fixing L, condition (i) above 
requires a{r) to be less than P + k elements in {a{ai), . . . , o"(a„_|_/3)}. Conditions (i) and 
(ii) together force exactly (3 + k + 771 elements of {l,...,m + n + /3 + l} to be more than 
(j(r). Thus the number of bijections a forming a valid pair {L,a) with this fixed L is 
{n + m + (3 + 1)1/ {k + m + j3 + 1) and the total contribution of all pairs (L, a) for this 
fixed L is (— l)"~'^(n + m + P + 1)1/ {k + m + P + 1). Let us call the pairs (L, a) satisfying 
(i) and (ii) valid pairs. 

Let us now eliminate the contribution of some valid pairs (L, a) by introducing the 
following involutions for i = 1,2, . . . ,n + p. If the Delannoy path L contains a northeast 
step from some {i — 1, y) to [i, y + 1) and a{r) < cr(aj) then we define Ti{L, a) as {L', a) 
where L' is obtained from L by replacing the northeast step from {i — 1, y) to {i, y + 1) with 
an east step from {i — 1, y) to {i, y) followed by a north step from {i, y) to {i, y+ 1). (Note 
that the condition cr(r) < a{ai) guarantees that [L', a) is a valid pair.) If L contains an 
east step from (z — 1, y) to (z, y) followed by a north step from {i, y) to (z, y+ 1), we define 
ri(L, 0") as {L', a) where L' is obtained from L by replacing the sequence of one east step 
and one north step from (i, y) to {i,y + l) by a diagonal step from {i, y) to (z, ?/ + 1). (Note 
that condition (i) for {L,a) requires cr(r) < a{ai) in this case.) In all other situations 
we define Ti{L,a) := {L,a). Obviously, the r^'s are involutions on the set of valid pairs, 
and they commute pairwise, since they change disjoint parts of the underlying Delannoy 
paths. Thus they induce a Z2'''^-action on the set of valid pairs. If Ti{L, a) 7^ {L, a) then 
the weight of Ti{L, a) is the negative of the weight of (L, a). Thus the weight of all (L, a) 
pairs that belong to the same Z2^^-orbit cancels, unless the orbit in question contains a 
single fixed point. 

We have shown that 



Jo 

equals the total weight of all valid pairs {L,a) satisfying Ti{L,a) = {L,a) for i = 
1,2, . . . ,n + p. These pairs may be characterized by the following two conditions: 

(iii) cr(r) > cr(aj) holds for all i such that there is a northeast east step in L from 
{i — l,y) to {i, y + 1) for some y; 

(iv) there is no east step immediately followed by a north step. 

Condition (iv) is a statement about the Delannoy path L only. It may be stated equiva- 
lently by requiring that the only way to go before the first, after the last, or between two 
consecutive northeast steps is to use all the north steps before all the east steps. Such 
a Delannoy path is uniquely determined by the set of its diagonal steps. A sequence of 
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diagonal steps can be completed to a Delannoy path satisfying (iv) if and only if the 
first coordinates and the second coordinates of the starting points both form a strictly 
increasing sequence. Thus there are Delannoy paths satisfying (iv). Given a 

Delannoy path L satisfying (iv) having k diagonal steps, conditions (i), (ii), and (iii) are 
equivalent to stating that a{r) = k + 1, a{ai) < k for all i such that there is a northeast 
east step in L from [i — 1, y) to (z, y + 1) for some y + 1, and the image of all remaining 
elements under a belongs to {k + 1, k + 2, . . . ,n + m + /? + l}. There are k\{n + m + f3 — k)\ 
ways to find such a a. 

Therefore we obtain the following equality: 
(19) {n + m + p + l)\- j\"'+^-Pj^^'''\x) dx = J2(-'^)''(^l'^^ (l) ■ k\{n + m + P - k)\ 



k / \k 



We are left to show that the right hand side is zero for m < n. This is known in the theory 
of rook polynomials (see Section H]), so here we indicate a short proof for completeness 
sake only. The right hand side is (m + /?)! times 

p(m) := X^(-l)'(^) (n+/3-A;+l),(m+/5)„_fc = (-1)"^^ Q {n+(3-k+l)u{-m-(3-n+k)n. 

k=0 ^ fc=0 ^ ^ 

Consider p{m) as a polynomial function of m. The number (— l)"p(— m) is then the 
number of ways to select a /c-element subset of an n-element set and injectively color its 
elements using n + (5 colors, then color the remaining n — k elements injectively, using a 
disjoint set of m — /5 — 1 colors. Thus we have 

f n + m — l\ fn — m — 1 

(—1) p[—m) = implying p[m) = (—1) 

\ n J \ n 

This is zero for m < n. □ 

The rest of the proof of Theorem 13.21 is a consequence of the following Proposition. 
Proposition 3.4. For all a,/? G the shifted Jacobi polynomials satisfy 



[X - l)"Pi°'^)(x) = f2 (") (-l)'a;"-'Pf '"+^-')(x). 



Proof. First we show using Proposition 12.61 that (x-l)°Pi"'^^(x) is the total weight of 
all Delannoy paths from (0, 0) to {n + a + P,n + a) subject to the restriction that none 
of the first a steps is an east step. Here each east step contributes a factor of 1, each 
north step contributes a factor of x and each northeast step contributes a factor of —1. 
Indeed, if such a Delannoy path has k north steps then it has n + a — k northeast steps 
and P + k east steps. There are ("'^") ways to determine the order of the north and the 
northeast steps among themselves and then there are ("^^) ways to determine the order 
of the P + k east steps with respect to the n other steps (we subtracted a from n + a for 
the first a steps which can not be east steps). 
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Every Delannoy path subject to the above restriction may be uniquely decomposed as 
follows. The first a steps form a Delannoy path with no east step from (0, 0) to some 
(i, a) where < z < a is the number of northeast steps among the first a steps. These 
contribute a factor of (— Given z, there are (") ways to determine of the order 
of the first a steps. The rest of the path is an unrestricted Delannoy path from (z, a) to 
(n + a + /?, n + a), shifting its start to the origin yields a Delannoy path from (0, 0) to 
(n + a + /? — i, n) . This part of the path may be chosen independently among all Delannoy 
paths from (0,0) to (n + a + /5 — z, n) and, by Corollary 13.11 their total contribution is 

Using Proposition 13.41 we may conclude the proof of Theorem 13.21 as follows. Assume 
a, /3 G N and m < n. Then 

[' x"" ■ Pi'^'^^x) ■ (1 - dx = (-1)° [' x'^^^ ■ V f (-l)^x"-^P^°'°+^-'Hx) dx 







1=0 

All integrals in the last sum are zero by Proposition 13. 3[ 

4. Connections to the classical theory of rook polynomials 

For /5 = equation f|T9|) takes the form 

(20) (n + m+1)!- / x"* ■ P„(a;) = V(-l)'= W ■ k\{n + m - k)\ 

The right hand side here is £(x'"/„(x)) where /„(x) is the Laguerre polynomial defined in 
(fT2|) . The orthogonality of these Laguerre polynomials with respect to the inner product 
induced by £ is equivalent to stating 

/•oo 

Cix^'lnix)) = / x™Z„(x)e-^ dx = for m < n. 
Jo 

Thus the combinatorial orthogonality of the shifted Legendre polynomials reduces to 
the combinatorial orthogonality of the Laguerre polynomials used in the theory of rook 
polynomials. At a hasty look, the connection is "almost obvious" visually, since in both 
combinatorial situations the term (2) expresses the number of ways to select k rows and 
k columns on some n x n "chess-board", as the rows and columns where diagonal steps, 
respectively, rooks are located. The main difference is in the next step of the selection 
process: in the case of the valid pairs (L, a) that are fixed points of the involutions Tj 
the diagonal steps must follow a strictly increasing order in the picture, and a factor of 
k\ arises by selecting the restriction of an injective map a to the columns associated to 
the diagonal steps. In rook polynomial theory, however, the selected rooks may have any 
permutation pattern, thus a factor of k\ arises in a completely different way. It should be 
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also noted that the condition m < n was not used in the proof of (fTOjl . thus we may state 
in general 

(21) (n + m + 1)! ■ / x"" ■ Pn{x) dx = x'"/„(a;)e-^ dx for all m,nEN. 

Jo Jo 

For general /3 G N, equation (fT9l) may be rewritten as 

(22) (n + m + /3 + 1) ! ■ / ■ P^^'^) (x) c/x = / x™/(f ^ (x)a:^e-^ c^x for all m, n G N. 
Here 

n 

/f(x) :=5:(-l)'= 
fe=o 

is the n-th generalized Laguerre polynomial associaied to the rectangular board x [n] 

(as a subset of itself). These generalized Laguerre polynomials are related to the usually 
normalized generalized Laguerre polynomials L^f^(x) by the formula 

(x) = (-l)"n!L(f)(x), 

and form an orthogonal basis with respect to the weight function x^e~^. 

Remark 4.1. The reader familiar with the work of Foata and Zeilberger [10] should also 
note that in their work too, showing the orthogonality of the Jacobi polynomials is in a way 
reduced to the study of a combinatorial model for the Laguerre polynomials (presented 
in [2]). The question naturally arises whether a closer relation could be found between the 
orthogonality proof of Foata and Zeilberger [10] and the weighted lattice path enumeration 
proof presented in this paper. Finding such a connection depends on surmounting the 
following difficulty: in the work of Foata and Zeilberger the parameters a and /5 are weights 
(associated to bipermutations) and left as variables, whereas in the present model a and 
P need to be a natural numbers, since the coordinates of the endpoints of the lattice paths 
considered are expressed in terms of these parameters. 




5. Shifted Jacobi polynomials with zero first parameter and negative 

integer second parameter 

Using weighted Delannoy numbers it is easy to connect most shifted Jacobi polynomials 
Pn'^'^^ (x) with negative integer (3 to the Jacobi polynomials with positive integer (3. Indeed, 
we have the following formula. 

Proposition 5.1. For /5 G N and n > f3 we have 

Pi''-^\x)=x^h-p{x). 

Proof. By Corollary 13. II we have 

p(0,/3)/ N _ jl,x,-l 
n — "'n+8,n- 
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Dividing the weight of all north and northeast steps by x we may take out a factor of 
and get 

p{0,l3)( \ _ n,l,l-l/x 

Since each east and each north step has the same weight, we may swap the horizontal 
and vertical axis and get 

71 K-^J — "'n,n+l3 ■ 

Finally, multiplying both sides by x^ and using the arising factor of x"^'^ on the right 
hand side to multiply the weight of each north and each northeast step by x yields 



>(0,/3) 



l,x,-l 
+/3- 



Here the right hand side is equal to Pn^p {x) by Corollary 13.11 The statement follows by 



replacing n with n — (5. 



□ 



Proposition 15.11 does not apply to Pn' (x) when n < /3. Table [U contains these 



polynomials and 



(0,-6), 



x) for 13 = Q. 



n 





1 


2 


3 


4 


5 


6 


P^''-'\x) 


1 


5 - 4x 


3x2 _ + 10 


3x2 _ + 10 


5 - 4x 


1 


x^ 



Table 1. P 



(0,-6), 



x) for n < 6 



The list of polynomials exhibits a certain symmetry which is most easily shown for the 
transformed Jacobi polynomials P^'^\2x + 1). By (fTOj) . for the values of n allowed by 
Theorem 11.11 these polynomials are equal to the Romanovski polynomials b!^'^\x). 

Definition 5.2. We extend the definition of the Romanovski polynomials R^'^\x) to all 
values of n using ( fTOil . 



As an immediate consequence of (jlj) and ( |TOl) we obtain the following. 
(23) 



n + a + (3 + + a\^^^ 



i=o ^ J / \ J 

In particular, substituting a = and —(3 for (3 yields 



i?f-^)(x) 
which may be rewritten as 
(24) R^n'~^\x) = 



E 

j=0 



(3 — 1 — n\ fn 
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Corollary 5.3. For P eN and < n < P — 1 we have 
implying also 

Pi''-'\^) = P'^r^iix) and p(o,-«(^) = p(°j-f) (x). 

Because of Corollary 15. 3[ the Jacobi polynomials {Pn~^\x)}n>o or, equivalently, the 
polynomials {Rn' ^\x)}n>Q can not form an infinite orthogonal polynomial sequence. 
There is a chance to having such a sequence only up to n = — 1)/2J. Theorem 11.11 
guarantees having a finite sequence of orthogonal polynomials up to almost that degree. 
Indeed, substituting a = and replacing (3 with —(3 in Theorem 11.11 the inequality 
a + j3 + N + 1 <0 becomes N < (3 — 1 which, for integers, is equivalent io N < (3 — 2. 
Thus, Theorem 11.11 states the orthogonality of the polynomials Rn~^\x) of degree at 
most {(3-2)/2 only. 

Corollary 5.4 (Askey-Romanovski). For any positive integer (3 > 2, the polynomials 
(x) : n < [(/? — 2)/2j} form a finite sequence of orthogonal polynomials. 



We will provide a combinatorial proof of a strengthened version of Corollary 15.41 that, 
for odd values of /5, allows adding the polynomial -R(°'_f)/2(^) to finite sequence of or- 
thogonal polynomials, thus filling the only remaining gap. Note that the linear functional 

POO 

(25) Cp{x^) := / x''{l + x)"^ dx for < A; < /3 - 2, 

Jo 

implicitly used in Theorem ll.il is by definition P, 6.2.1] an evaluation of the beta function 
B{z,w). Thus we obtain 

(26) Cfsix") = B{k + l,(3-l-k) = j'^ ij^t) ~ ^^"'^ 

Since B{z,w) = r{z)r{w)/r{z + w) (see [H 6.2.2]) and r{z + 1) = z\ holds for all z G N 
(see by [H 6.1.15]), we also have 

(27) = ^nl^Tyr^ forO<fc</5-2. 

Keeping these formulas in mind, we are ready to present a combinatorial proof of our 
extended version of Corollary 15. 4[ 

Theorem 5.5. Let (3 > 2 he any positive integer and let Cp be the linear functional 
defined on the vector space {p{x) G C[x] : deg(p) < /3 — 2} by fl25l) . Then the Romanovski 
polynomials {Rn'~^\x) : < n < (/3 — 2)/2} form an orthogonal basis in the with respect 
to inner product {f,g) := Cp{f-g) on the vector space {p{x) G C[x] : deg(p) < (/5 — 2)/2}. 
For odd (3 we may extend Cp and the induced inner product to polynomials of degree 
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at most {(3 — l)/2 hy making Cf3{x^ ^) large enough to make the determinant of the 
{P + l)/2x {p + l)/2 matrix 

I Cp{x') Cp{x^) ■■■ Cp{x^P-^y^) \ 

positive. The polynomial Rf^^Iiy2i^) "^^?/ then be added to the orthogonal basis. 

Proof. The fact that the bihnear operation {f,g) := ■ g) in an inner product is a 

consequence of the integral form of Cj^i^x^) stated in (1261) . Using this equation it is easy 
deduce 

(28) (/, 9) = (y^) ■ 9 (y^) ■ (1 - tf~' dt if deg(/) + deg(^?) < (3 - 2. 

Note that the condition deg(/) + deg{g) < (3 — 2 imphes that the degree of (1 — t) in the 
denominator of f{t/{l — t)) ■ g{t/{l — t)) is at most /3 — 2, and so the integrand on the 
left hand side is a polynomial function. Thus the bilinear function is in deed an inner 
product. To prove the first half of the theorem, it is sufficient to show 



(x™, =0 for m < n < ^ 



-2 



2 



§-2 
2 ■ 



By ( 1271) this is equivalent to 

^ ~ ^ ~ (^^ {-iy{m + j)\{p-2-m- j)\ = for m < n < ^ 

Dividing both sides by (/3 — 1 — n)\m\{n — m — 1)\ yields the equivalent form 

The left hand side of f l29|) is the total weight of all triplets {X, A, B) where 

(i) X is a subset of {1,2,..., n}; 

(ii) A = {oi, . . . , Qm} is an m-element multiset such that each Oj belongs to X U {0}; 

(iii) B = {bi, . . . , bn-m-i} is an (n — m — l)-element multiset such that each bj belongs 
to {1,...,(3 -n}\X. 

The weight of the triplet (X, A, B) is (—1)'"^'. Indeed, for a fixed j there are (") ways to 
select X. Once X is fixed, there are 

\j + 1) + m — l\ / m+ j 
m J \ m 

ways to select the multiset A on X U {0} and 

\(3 — n — j) + {n — m — 1) — l\ //3 — m — 2 — j 
n — m — 1 J \ n — m — 1 
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ways to select the multiset i? on {1, . . . , /3 — n — 1} \ X . Let us now fix the multisets A 
and B first. Since 1^41 + \B\ = n — 1, there is at least one c G {1, . . . ,n} that does not 
appear in A, nor in B. Then, for a subset X C {1, . . . ,n}\ {c}, the triplet {X,A,B) 
satisfies the above conditions if and only if the triplet {X U {c}, A, B) does. The weight 
of these triplets cancel each other. 

For the case when (3 is odd, observe first that the above proof of orthogonality can be 
adapted to m < n < (/5 — l)/2 without any changes since we only need to use the value of 
jCj3{x^) for k < (3 — 2. The only thing left to show is that we can set the value of £/3(a;^~^) 
in such a way that we get a positive definite inner product. It is well known that a bilinear 
form defined by a symmetric real matrix is positive definite if the principal minors of the 
matrix are positive. Thus we only need to make sure that the principal minors of Mp are 
positive. The fact that the restriction oi Cp io the vector space of polynomials of degree 
at most {(3 — 3)/2 is positive definite implies the positivity of all principal minors of Mg 
except for det(Mg). Finally, consider Cpi^x^) as a variable, and expand det(Ma) by its 
last column. The determinant becomes a linear function of whose coefficient is the 

(/3 — 3)/2 X {13 — 3)/2 principal minor of Afg, a positive number. Therefore det(Mg) > 
holds for any sufficiently large Cj3{x^). □ 

6. Weighted Schroder numbers and iterated antiderivatives of Legendre 



We define a Schroder path from (0, 0) to (n, n) as a Delannoy path not going above the 
line y = X. In analogy to Definition 12.11 we may introduce weighted Schroder numbers as 
follows. 

Definition 6.1. Let u, v, w be commuting variables. We define the weighted Schroder 
numbers s]^'""'^ as the total weight of all Schroder paths from (0,0) to {n,n), where each 
east step (0, 1) has weight u, each north step has weight v, and each northeast step has 
weight w. The weight of a lattice path is the product of the weights of its steps. 

Introducing the Schroder polynomials by 
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Sn{x) : 




for n > 0, 



Proposition 6.2. For n > 1 the Schroder polynomials are given by 



{n + l)x 




Proof. First we show that the Schroder polynomials satisfy 



(31) 




for n> 1. 
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Let j G {0, . . . , n} be the number of east steps in a Schroder path, the number of norths 
steps must be the same, the number of northeast steps is {n — j). The path is a Schroder 
path if and only if at any given step the number of north steps thus far does not exceed 
the number of east steps. The number of ways to arrange j east steps and j north steps 
satisfying this criterion is the Catalan number /(j + 1), and there are (^I'i^) ways to 
determine the order of the northeast steps with respect to the east and north steps. Since 

'2j\/n + j\^/ n+j \^^n + j\/n^ 
j J\n- jj \J, J, n-jj \ n J\j^ 
1 / n\ 1 f n + 1 



and 



j + 1 \j J n + l\j + 1 



we may rewrite (13T|) as 



On the other hand. Proposition 12.61 gives 



for n > 1. 



n+l 



[X 



k=0 



l)p(W)^3,^ = ^(_l^n+l-fc^fc -yM - forn > 1. 



n 

Observe that for n > 1 the constant term in the last sum is zero, thus we may write 

'n\ fn — 1 + k^ 
n 



n+l 
k=l 



The statement now follows from substituting j ' + 1 into k in the last equation and com- 
paring the result with (132|) . □ 



In analogy to Lemma 12.31 we have the following 

Lemma 6.3. The weighted Schroder numbers are linked to the Schroder polynomials by 
the formula 



w 



J 



As a consequence of Proposition 16.21 and Lemma 16.31 we obtain 

Corollary 6.4. Forn > 2, the weighted Schroder numbers are linked to the shifted Jacobi 
polynomials by the formula 



n+lV uv/ \ w / 



n 



In particular, substituting u = v = w = 1 gives that the ordinary Schroder numbers 
are given by 
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Using the swapping rule ([5]) we obtain the following remarkable variant of ([2]). 

(33) s^ = _J_.pM,i)(3) forn>l. 

The same swapping rule for shifted Jacobi polynomials (ITSjl allows to rewrite Corollary 16 .41 
as 

(34) = ^ f 1 + -) f^ + i) forn > 1. 
^ ^ " n + 1 V uvJ \w J ~ 

Central Delannoy numbers and Schroder numbers satisfy an obvious recursion formula, 
which may be generalized to their weighted variants as follows. 

Proposition 6.5. The weighted central Delannoy numbers and the weighted Schroder 
numbers satisfy the following formula. 



n-1 



f]u,v,w _ ^ m,v,w u,v,w , ju,v,w 



k=0 

Proof. For any Delannoy path from (0,0) to (n, n) there is a unique {k,k) on the path 
where the path contains a lattice point of the form (i,i) for the last time before {n,n). 
Consider the contribution of all Delannoy paths associated to the same fixed (fc, k), where 
k G {0, 1, . . . , ra — 1}. The part of all such paths up to (A;, k) contributes a factor of d^'l'^ . 
The rest must begin with an east step and end with a north step, contributing a factor 
of uv. For k < n — 1 the remaining steps contribute a factor of 2s^'^^j^. The factor of 2 
represents the choice of staying strictly above or below the line y = x for the rest of the 
path. For k = n — 1, besides the contribution of 2uvSq'^''^ = 2uv the possibility of adding 
a northeast step from (n — l,n — 1) to (n, n) arises: this possibility is accounted for by 
the term □ 

Using Corollary 13.11 and (|5U]) we obtain: 
Corollary 6.6. The shifted Legendre polynomials satisfy the recursion formula 

n-1 



Pn{x) = 2x^Pfc(x)5'„„fe_i(x) - Pn-l{x) 



k=0 

Using Proposition 16.21 the last corollary may be rewritten as follows. 
Corollary 6.7. We have 

n-2 ^ 

Pn{x) = 2 V P,{x)^^Pi':l^\ix) + {2x - l)P„_i(x) and 

k=0 

n-2 

Pn{x) = y^Pk{x)^^Pi'rk!.\ix) + xPn-lix) forn> 1. 

Tl rC 

k=0 
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Comparing (!3T]) with ([8]) we find immediately that the Schroder polynomials are closely 
related to the the antiderivatives of the shifted Legendre polynomials: 

1 r ~ 

(35) Sn(x) = - / Pn-i(t) dt holds for n > 1. 

X Jo 

By Proposition 16.21 this implies 

(36) (x - l)P^^'-^\x) = / Pn(t) dt for n > 1. 

n+l Jo 

This equation may be generalized to the following statement. 

Proposition 6.8. Let n and a he positive integers. Applying the antiderivative operator 

f{x)^ r fit) dt 







to Pn{x) exactly a times yields the polynomial {x — l)"Pi"' °'\x). 

Proof. By Proposition 12.61 we have 



Here the binomial coefficient ) is zero for k < a. Thus, after introducing j := k — 



a, 



we may write 



X - l)°P^"'-")(a;) = n^)' 



n \ ( n — ] 



3=0 

Thus 

Taking the derivative on both sides yields a right hand side of the same form, with the 
value of a decreased by one. Thus we obtain 

d {x-\YPt-''\x) (x-l)"-ipi"~''"("-'))(x) ^ 

r = r tor a > 1. 

dx [n + 1)q, [n + 1)q_i 

The statement follows from applying this observation a times and from the fact that 
P^°'-°)(0) = holds for n, a > 0. □ 

Proposition 16.81 has an interesting consequence for the Narayana polynomials Nn{x) 
defined by 



xK 
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As it was shown by Mansour and Sun [23l Theorem 2.1], the Narayana polynomials are 
related to the antiderivatives of the shifted Legendre polynomials via the formula 

(39) A^„,(x) = {x- 1)"+^ / P„(2t -l)dt={x- 1)"+^ Pn{t) dt. 

Jo Jo 

Direct substitution of Proposition 16.81 into (1391) yields 



n + lVx — 1 / Vx — 1 



Corollary 6.9. The Narayana polynomials Nn{x) are connected to the shifted Jacobi 
polynomials P^^'~^\x) by the formula 



n+1 \x — 1 J 



1 am indebted to Yidong Sun for helping me realize this connection. 

We conclude this section by observing that the polynomials {5'„(x)}„>o almost form an 
orthogonal polynomial sequence. In fact, it is possible to show that the monic polynomials 

, , 1 x-1 



Pn[X) 



satisfy Favard's recursion formula (ITTl) 



Unfortunately, substituting n = 2 yields A2 = thus the moment functional will not be 
quasi-definite. 



7. Concluding remarks 



One of the main observations in this paper is that the connection ([T]) between the central 
Delannoy numbers and Legendre polynomials may be extended to all Delannoy numbers 
and the Jacobi polynomials {Pn°''^\x)}n>o via ([2]). An eerily similar relation was found 
in fTEl (3.2)] between the modified Delannoy numbers dm,n and the Jacobi polynomials 
{P^''^\x)}n>o- Formula (3.2) in [15] may be restated as 

(40) rf„+^,„ = Pf^)(3) for/3GN. 

The modified Delannoy number dm,n is the number of lattice paths from (0, 0) to (m, n+1) 
whose steps belong to N x P. (Here P denotes the set of positive integers.) It would be 
worth exploring whether the weight enumeration of modified Delannoy paths (defined as 
having steps from N x P) could also be done using shifted Jacobi polynomials, and whether 
there is a duality between the two theories that extends the ( "signless" ) swapping of the 
parameters observed above in ([2]) and (140|) . 
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The fact that the (shifted) Legendre and Jacobi polynomials are orthogonal with respect 
to some weight function is also a consequence of Favard's theorem. The question naturally 
arises whether Viennot's combinatorial theory [3lj could be applied to find a moment 
functional with respect to which they are orthogonal. Laguerre polynomials (their usual 
monic variant) have a very nice combinatorial interpretation in Viennot's work. For 
(the monic version of the) Legendre and Jacobi polynomials, however, the coefficients 
{cn}n>i and {A„}„>i are not all integers. For example, the monic variant of the Legendre 
polynomials is given by 

PnyX) . ^2n-j 

Favard's recursion formula (fTTD takes the form 



Pn{x) = XPn-lX - — — —Pn-2{X). 

[2n — lj[2n — S) 

Thus some weights — c„ and — A„ used on Viennot's Favard paths have to be fractions, 
and the formula expressing the moments as total weights of Motzkin paths seems very 
hard to evaluate. For Legendre polynomials, the horizontal steps will have zero weight, 
the northeast steps (1, 1) will have weight 1, the southeast steps (1,-1) will have weight 
A;^/(4fc^ — 1) if they start at a point whose second coordinate k. Using the fact that 
Legendre polynomials are orthogonal with respect to the inner product 

{f^a) ■= I f{x) ■ g{x) dx, 



J -I 

we know that the moments must be a constant multiple of 

1 , f if X is odd, 
— dx = < 2 -t ■ 

Experimental evaluation of Viennot's Motzkin paths for the first few even values of n 
shows that the total weight of these paths is l/(n + 1). A direct combinatorial proof of 
this fact would be desirable. It should also be noted that the non-monic variant 

2"(2n-l)!!P^(x) 

of the Legendre polynomials satisfies the recursion formula 

(41) qn{x) = {2n - l)xg„_i(x) - {n - l)^g„_2(x) 

with integer connecting coefficients. Perhaps a non-monic version of Favard's theorem is 
worth considering, to which Viennot's combinatorial proof could be extended by replacing 
the weight x on the short vertical steps in his Favard paths by a weight of {2n — l)x where 
n is determined by the height of the step. 

We should also mention the combinatorial model of the Jacobi polynomials constructed 
by Foata and Leroux [8], whose study was continued by Leroux and Strehl [21], [30], [31] . 
This model was never used to prove the orthogonality of the Jacobi polynomials, and it 
seems hard to be used for such a purpose. On the other hand, finding an interpretation 
for the Delannoy numbers in it seems to be an interesting question for future research. 
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As seen in Section HI our current work relates the Legendre and Jacobi polynomials to 
the classical theory of rook polynomials. The level of this connection is not satisfying, since 
we only found a binomial identity to which the orthogonality of the Laguerre polynomials 
(rook theory variant) and the orthogonality of certain Jacobi polynomials may both be 
reduced. It would be desirable to find results directly relating colored Delannoy paths 
and rook placements as combinatorial structures. 

Our proof of the orthogonality of the Jacobi polynomials P^'^\x) for nonzero a de- 
pends on the algebraic formula given in Proposition 13.41 A more direct combinatorial 
extension to the case a > would be desirable, which can perhaps be found once the "dual- 
ity" between Pn^'^\x) and Pn'^\x) is better understood. Although Proposition 16.81 seems 
to give a "philosophical reason" to focus on the polynomials {x — 1)°' P^~°'\x) / {n + 1)q,, 
these have a combinatorial interpretation for a = 1 only, for higher values we start hav- 
ing non-integer coefficients. It is remarkable though that the same polynomials [x — 
l)<^Pt""\x), without being divided by {n + l)a play a key role in our work and have a 
combinatorial interpretation for all a G N. 

Finally we wish to note that there is a yet to be explored connection with the theory 
of Riordan arrays (see Sprugnoli [28] for definition and uses). The weighted Delannoy 
number d^^^ is the coefficient of in [u + wt)"^/{l — vt)'^'^^. An immediate conse- 
quence of this observation is that the n-th row fc-th column entry in the Riordan array 
(1/(1 — vt),t{u + wt)/{l — vt)) is d^'^^i^. The numbers dj^^^^ appear as entry A1016195 
in Sloane listing the entries of the Riordan array (1/(1 — 2t), t(l — t)/(l — 2t)). Our 
results should allow to write summation formulas for Jacobi polynomials using the theory 
of Riordan arrays. 
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